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Abstract For the case of Tycho’s supernova remnant (SNR) 
we present the relation between the blast wave and contact 
discontinuity radii calculated within the nonlinear kinetic 
theory of cosmic ray (CR) acceleration in SNRs. It is demon¬ 
strated that these radii are confirmed by recently published 
Chandra measurements which show that the observed con¬ 
tact discontinuity radius is so close to the shock radius that it 
can only be explained by efficient CR acceleration which in 
turn makes the medium more compressible. Together with 
the recently determined new value E sn — 1.2 x 10 51 erg of 
the SN explosion energy this also confirms our previous con¬ 
clusion that a TeV y-ray flux of (2 — 5) x I 0“ 13 erg/(cm 2 s) is 
to be expected from Tycho’s SNR. Chandra measurements 
and the HEGRA upper limit of the TeV y-ray flux together 
limit the source distance d to 3.3 < d < 4 kpc. 

Keywords (ISM:)cosmic rays - acceleration of particles 
- shock waves - supernovae individual(Tycho’s SNR) - 
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1 Introduction 


to individual SNRs f9lH0ll 11 11281 has demonstrated its power 
in explaining the observed SNR properties and in predicting 
new effects like the extent of magnetic field amplification, 
leading to the concentration of the highest-energy electrons 
in a very thin shell just behind the shock. 

Recent observations with the Chandra and XMM-Newton 
X-ray telescopes in space have confirmed earlier detections 
of nonthermal continuum emission in hard X-rays from young 
shell-type SNRs. With Chandra it became even possible to 
resolve spatial scales down to the arcsec extension of indi¬ 
vidual dynamical structures like shocks j27ll22l [3l. The fila¬ 
mentary hard X-ray structures are the result of strong syn¬ 
chrotron losses of the emitting multi-TeV electrons in am¬ 
plified magnetic fields downstream of the outer accelerating 
SNR shock 1 1271110111211301 . Such observational results gain 
their qualitative significance through the fact that these ef¬ 
fective magnetic fields and morphologies turned out to be 
exactly the same as predicted theoretically from acceleration 
theory. 

This theory has been applied in detail to Tycho’s SNR, 
in order to compare results with the existing data l28ll30l . 
We have used a stellar ejecta mass M e j — 1.4M 0 , distance 
d — 2.3 kpc, and interstellar medium (ISM) number den¬ 
sity i\'h — 0.5 H-atoms curl For these parameters a to¬ 
tal hydrodynamic explosion energy E sn — 0.27 x 10 51 erg 
was derived to fit the observed size R s and expansion speed 
V s . A rather high downstream magnetic field strength B,[ « 
300 Li G and a proton injection rate rj — 3 x 10~ 4 are needed 
to reproduce the observed steep and concave radio spectrum 
and to ensure a smooth cutoff of the synchrotron emission 
in the X-ray region. We believe that the required strength of 
the magnetic field, that is significantly higher than the MHD 
compression of a 5 fiG ISM field, has to be attributed to non¬ 
linear field amplification at the SN shock by CR acceleration 
itself. According to plasma physical considerations Il24l l5ll4l. 
the existing ISM magnetic field can indeed be significantly 
amplified at a strong shock by CR streaming instabilities. 

After adjustment of the predictions of the nonlinear sphe- 
rically-symmetric model by a physically necessary renor¬ 
malization of the number of accelerated CR nuclei to take 
account of the quasi-perpendicular shock directions in a SNR, 


Cosmic rays (CRs) are widely accepted to be produced in 
SNRs by the diffusive shock acceleration process at the outer 
blast wave (see e.g. IU71U3l l6l i2Tl|23l for reviews). Kinetic 
nonlinear theory of diffusive CR acceleration in SNRs |[7j[8] 
couples the gas dynamics of the explosion with the particle 
acceleration. Therefore in a spherically symmetric approach 
it is able to predict the evolution of gas density, pressure, 
mass velocity, as well as the positions of the forward shock 
and the contact discontinuity, together with the energy spec¬ 
trum and the spatial distribution of CR nuclei and electrons 
at any given evolutionary epoch t, including the properties 
of the nonthermal radiation. The application of this theory 
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very good consistency with the existing observational data 
was achieved. 

Using Chandra X-ray observations 1321 have recently 
estimated the ratio between the radius R c of the contact dis¬ 
continuity (CD), separating the swept-up ISM and the ejecta 
material, and the radius R s of the forward shock. The large 
mean value R c /R s — 0.93 of this ratio was interpreted as evi¬ 
dence for efficient CR acceleration, which makes the medium 
between those two discontinuities more compressible. 

Here we present the calculations of the mean ratio Rc/Rs, 
which are the unchanged part of our earlier considerations 
II281I30II . and demonstrate that these results, which are in 
fact predictions, fit the above measurements very well. Since 
our calculations have been made in spherical symmetry they 
concern a priori an azimuthally averaged ratio R c /R s - We 
shall extend them by taking the effects of the Rayleigh-Taylor 
(R-T) instability of the CD into account. We shall in ad¬ 
dition discuss a physical mechanism that leads to the ob¬ 
served azimuthal variations of R c /R s . Finally we shall take 
the recent determination of the mechanical energy output 
E sn ~ 1.2 x 10 51 J2) which results from the theory of de¬ 
layed detonations in the physics of type la SN explosions 
1191201 , in order to predict a range for the y -ray spectrum 
that is consistent with the existing upper limits of high en¬ 
ergy y -ray fluxes from Tycho's SNR. 


2 Results and Discussion 

Fig. 1 and partly Fig.2 show the calculations of shock and 
CD related quantities which were part of our earlier consid¬ 
erations M28II30I . The calculated shock as well as CD radii 
and speeds are shown as a function of time for the two dif¬ 
ferent cases of interior magnetic field strengths B,\ — 240 liG 
and B,\ — 360 pG considered, together with the azimuthally 
averaged experimental data available at the time. 

According to Fig.la Tycho is nearing the adiabatic phase. 
To fit the spectral shape of the observed radio emission we 
assumed a proton injection rate 7] = 3 x 10 -4 . This leads to 
a significant nonlinear modification of the shock at the cur¬ 
rent age of t — 428 yrs. A larger magnetic field lowers the 
Alfvenic Mach number and therefore leads to a decrease of 
the shock compression ratio, as seen in Fig. lb. The result is a 
total compression ratio a — 5.7 and a subshock compression 
ratio o s — 3.5 for B& — 240 pG. In turn o — 5.2, o s = 3.6, 
for Bi — 360 pG. 

Therefore, as can be seen from Fig.2, including CR ac¬ 
celeration at the outer blast wave, the calculated value of 
the ratio R c /R s for B,\ — 360 pG is slightly lower than for 
Rd = 240 pG. At the current epoch we have R c /R s ~ 0.90 
which is lower than the value R c /R s — 0.93 inferred from 
the observations. Qualitatively our result goes in the same 
direction as calculations by fPfll which modeled SNRs with 
a uniform specific heat ratio y e g < 5/3, for the circumstellar 
medium and the ejecta material alike. 

Projecting a highly structured shell onto the plane of the 
sky tends to favor protruding parts of the shell. Therefore the 
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Fig. 1 (a) Shock radius R s , contact discontinuity radius R c , shock 
speed V s , and contact discontinuity speed V c , for Tycho‘s SNR as func¬ 
tions of time, including particle acceleration; (b) total shock (cr) and 
subshock (<j s ) compression ratios. The dotted vertical line marks the 
current epoch. The solid and dashed lines correspond to the internal 
magnetic field strength B A = 240 pG and (B d = 360 pG), respectively. 
The observed mean size and speed of the shock, as determined by radio 
measurements Eg], are shown as well. 


average radius measured in projection is an overestimate of 
the true average radius. Analysing the amount of bias from 
the projection for the shock and CD radii 1(321 found a cor¬ 
rected "true” value R c /R s — 0.93 which is lower than their 
measured "projected average” value R c /Rs — 0.96, as a re¬ 
sult of the above geometrical effect. 

In turn, starting from a spherically symmetric calculation 
of the CD radius, as we do, one has to take into account that 
the actual CD is subject to the R-T instability. In the non¬ 
linear regime it leads to effective mixing of the ejecta and 
swept-up ISM material with "fingers” of the ejecta on top of 
this mixing region, which extend farther into the shocked gas 
than the radius R c predicted when assuming spherical sym¬ 
metry e.g. nMMHum Therefore our ratio R c /R s = 0.90, 
calculated within the spherically symmetric approach, has 
to be corrected for this effect in order to compare it with the 
measured value Rc/Rs — 0.93. In the case when all the fin- 
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Fig. 2 The ratio R c /R s of the radii of the contact discontinuity and 
the forward shock as a function of time. Solid and dashed lines corre¬ 
spond to the same two cases as in Fig. 1. The lines of all other styles 
correspond to the SN explosion energy E sn = 1.2 x 10 51 erg and four 
different distances (see the caption of Fig.3) Thin lines represent the 
values calculated in the spherically symmetric model, whereas the thick 
lines show the values R',./R s which contain the correction for the effect 
prod uced by the R-T instability. The experimental point is taken from 


gers have length l and occupy half of the CD surface, one 
would have a mean CD size R' c « R c + 0.5 1 which has to 
be compared with 0.93R s . According to the numerical mod¬ 
elling of ||3T1 , albeit without particle acceleration, the R-T 
instability allows fingers of ejecta to protrude beyond the 
spherically symmetric CD radius by 10%. The longest fin¬ 
gers of size l « 0.1 R, occupy less than 50% of the CD sur¬ 
face. However, in projection they stick out of the mixing re¬ 
gion, whose thickness is roughly 0.5/. This leads to a rough 
estimate of the corrected CD radius R' c — 1 .05 R, which has 
to be compared with the experimentally estimated value. 

The comparison of the corrected values R'/R s , accord¬ 
ing to our earlier calculations as well as for different assump¬ 
tions (see below) about the explosion energy and source dis¬ 
tance, with this experimentally estimated value R c /R s = 0.93 
(in Fig.2 we present that value with 2% uncertainties, ac¬ 
cording to [[32] ) shows quite good agreement (see Figs.2 and 
3) even if one takes into account some uncertainty in the 
quantitative determination of our correction factor R' c /R c , 
which in our view lies in the range 1.03-1.07. 

Another interesting peculiarity of Tycho’s shock struc¬ 
ture, which we would like to discuss here, is the quite ir¬ 
regular behaviour of the radius of the forward shock around 
the edge of the visible SNR disk, that is clearly seen in FigOJ 
Large shock distortions of this kind are not expected to result 
from the R-T instability. At first sight such a variation of the 
values of R s and R c can be easily attributed to fluctuations 
of the ambient ISM density and/or to an inhomogeneously 
distributed density and velocity field of the ejected matter. 
The local shock part, which encounters a lower ISM density 
or has a faster ejecta portion behind, propagates faster com¬ 
pared with the neighbouring shock pieces. This will lead to 
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Fig. 3 Top panel: the forward shock radius R s and bottom panel: the 
ratio R c /R s of the radii of the contact discontinuity and the forward 
shock as a function of azimutal angle (32) . The regions where R c /R s > 
0.99 are shadowed. 


the formation of local forward displacements of the CD and 
the forward shock whose number and relative sizes are de¬ 
termined by the specific structure of the ISM and/or ejecta. 
Since a larger shock compression ratio is expected for higher 
shock speed (see Fig. lb), one should also expect a smaller 
difference between R s and R c on the top region of each such 
displacement. This is exactly what is observed. 

Such arguments would give a good explanation for the 
observed picture only if CR injection/accelerationtook place 
uniformly across the entire shock surface. However, in our 
view the actual situation is expected to be more compli¬ 
cated - and physically more interesting. Efficient injection 
of suprathermal nuclear particles into the acceleration pro¬ 
cess takes place in those local shock regions, where the for¬ 
ward shock is quasi-parallel, and these regions are distributed 
over the shock surface according to the ambient ISM mag¬ 
netic field structure and occupy in total about 20% of the 
shock ll29l . As a result one would expect that only about 
20% of the local forward shock displacements efficiently 
accelerate nuclear CRs, and therefore display the extraor¬ 
dinary high ratio R c /R s . In reality the most extreme values 
Rc/Rs > 0.99 are observed on randomly distributed local 
shock regions, and the positions of these regions roughly co¬ 
incide with the positions of extended shock displacements. 
To explain such a correlation a strong physical connection 
between the shock speed and the efficiency of CR injec¬ 
tion/acceleration should exist: the local parts of the shock 
with large speed effectively produce CRs and vice versa. 

There are at least two physical processes which can re¬ 
solve the above problem. The first one gives high injection 
on the leading part of the displacements, if they are formed 
due to ISM and/or ejecta inhomogeneities and initially the 
injection was suppressed here as a result of a highly oblique 
magnetic field. Since after its formation the leading part of 
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the displacement develops a large curvature, it will have a 
significant portion which becomes quasi-parallel, and there¬ 
fore efficient CR injection/acceleration is expected over the 
whole leading part of such a displacement. 

The second factor can presumably itself lead to the for¬ 
mation of displacements, even in the case of a uniform ISM 
and a spherically symmetric ejecta distribution. A signifi¬ 
cant fraction of the internal energy behind a quasi-parallel 
part of the shock front is contained in CRs. Since the shock 
is expected to be modified, the CR spectrum is very hard and 
therefore the pressure is mainly carried by the CRs with the 
highest energies. With their high mobility these CRs can dif¬ 
fuse laterally into neighbouring downstream volumes, which 
are located behind shock surfaces not producing CRs. This 
diffusive loss leads to a decrease of the internal pressure 
behind the quasi-parallel parts of the shock. Therefore the 
corresponding downstream ejecta undergo less deceleration. 
As a result these ejecta accelerate relative to their surround¬ 
ings. This leads to the formation of a local displacement. 
One can expect that this also distorts the initially spherically 
symmetric shock. We therefore conclude that the effectively 
accelerating parts of the forward shock are among those re¬ 
gions, where the displacements occur, and they may actually 
delineate them. Predominantly fingers are situated on the ef¬ 
fectively accelerating parts of the shock surface. If this is 
correct it gives a unique experimental identification of the 
areas at the outer SNR shock which efficiently produce nu¬ 
clear CRs. 

Since the most outwardly displaced parts of the CD and 
the shock with a high ratio R c /R s have to be interpreted 
as the regions with efficient CR injection/acceleration, this 
opens the possibility to experimentally distinguish the shock 
areas with efficient CR production from those where CRs are 
not produced. A rough estimate shows that in Tycho's SNR 
the regions with extremely high ratios R c /R s >0.99 and the 
displacements occupy about 20% of the shock surface see 
Fig.4 of j32j. This corresponds to the theoretical expecta¬ 
tion. 

Our last point regards constraints which the recent re- 
evaluation of the mechanical output E sn — 1.2 x 10 51 erg to¬ 
gether with the HEGRA upper limit for the TeV y-ray flux 
m approximately impose on the distance and ambient den¬ 
sity for Tycho’s SNR. With this new £ sn -value we find a con¬ 
sistent fit for all existing data - the SNR size, its expansion 
rate, overall synchrotron spectrum and the filament struc¬ 
ture of the X-ray emission - like it was done for the previ¬ 
ously defined explosion energy E sn = 0.27 x 10 51 erg (see 
12811301 for details). In particular, the fit of R s and V s , which 
is of the same quality as in Fig.l for E sn — 0.27 x 10 51 erg, 
gives for each assumed distance d a rather definite value 
of the ISM number density Nn(d). A rough explanation of 
this numerical result is the following: since for a given ex¬ 
perimentally measured angular SNR size and its expansion 
rate the linear size R s and the speed V s scale proportion¬ 
ally to distance d, and since in the nearby (in time) Se¬ 
dov phase R s °c ( E sn /;\jj) 1 /<5 , the density Ah E sn /d 5 de¬ 
creases with increasing distance d. The hadronic y-ray flux 



Fig. 4 Spectral energy distribution of the y-ray emission from Tycho’s 
SNR, as a function of y-ray energy £y, for a mechanical SN explosion 
energy of E sn = 1.2 x 10 51 erg and four different distances d and cor¬ 
responding values of the ISM number densities Ah. All cases have a 
dominant hadronic compared to Inverse Compton y-ray flux. Experi¬ 
mental data are the upper limits of the HEGRA (H-CT; fil l and Whip¬ 
ple (W; fl5l l Cherenkov telescopes and the 95% confidence HEGRA 
A1ROBICC (HA; (25]) upper limit. 


Fy os A , s V s 2 A , f|/c/ 2 is then expected to scale as Fy ElJd\ 
The ejected mass is still assumed to be M c] — 1.4M 0 . We 
also find the same nuclear injection rate i)=3x 1(T 4 for all 
cases, and downstream magnetic field values B,\ « 400 /iG. 
At the same time, the linear size L of an X-ray filament in¬ 
creases proportional to d. Therefore the magnetic field strength 
B' d °c Z> 2 / 3 |fl2] , determined from the filament sizes (see In¬ 
troduction), decreases with d. 

In order to find the constraint on the distance d and the 
ISM density Ah, we then compare in Fig.3 the resulting y- 
ray spectral energy distribution with the HEGRA and Whip¬ 
ple upper limits at TeV energies. It is seen that all distances 
d <3.3 kpc are inconsistent with the HEGRA data. Dis¬ 
tances of d > 4 kpc are still consistent with the y-ray data, 
although there is a growing discrepancy between B^ and B' d : 
at d — 4.5 kpc B' d ~ 300 juG which is already consider¬ 
ably smaller than B& m 400 /iG. Therefore we believe that 
we can constrain the source distance also from above, d < 

4 kpc. In Fig.2 we also show the values for R c /R s and R' c /R s 
for the case E sn — 1.2 x 10 51 erg and these increased dis¬ 
tances. Within our approximate determination of R' c from R c 
they still agree with the Chandra data, in particular because 
the CR production rates are comparable. Our calculations of 
the y-ray emission lead us to predict that the new Northern 
Hemisphere TeV detectors should detect this source at TeV- 
energies in, predominantly, hadronic y-rays: the expected 
7T°-decay y-ray energy flux (2 — 5) x 10" 13 erg/(cm 2 s) ex¬ 
tends up to almost 100 TeV if the distance is indeed within 
the range 3.3-4 kpc. As a corollary the detection of a TeV 
signal is not only important by itself, but it is also crucial for 
the correct determination of all other key Supernova param¬ 
eters. 
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